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Abstract 

We extend the work of Fuchs, Painleve and Manin on a Calogero-like expres- 
sion of the sixth Painleve equation (the "Painleve-Calogero correspondence") to 
the other five Painleve equations. The Calogero side of the sixth Painleve equa- 
tion is known to be a non-autonomous version of the (rank one) elliptic model of 
Inozemtsev's extended Calogero systems. The fifth and fourth Painleve equations 
correspond to the hyperbolic and rational models in Inozemtsev's classification. 
Those corresponding to the third, second and first are not included therein. We 
further extend the correspondence to the higher rank models, and obtain a "multi- 
component" version of the Painleve equations. 
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I Introduction 



The so called Painleve equations are the following six equations discovered by Painleve 
and Gambier 0: 

. d 2 X 1 (\ 1 1 \ (d\V (\ 1 \ \d\ 

( p vi) -jw = 7T t + ^-r + ^— : Hi 



dt 2 2 VA A-l \-tJ\dtJ \t t-l X-tJdt 
A(A-l)(A-t) / /K 7 (t-l) , ft(t-l) 

H ^ TTo C* + TT7 + 7^ TTn + 



t 2 (t-l) 2 V A 2 (A-l) 2 (A-t) 2 
d 2 X ( 1 1 \ /dA\ 2 IdA 



^ dt 2 \2X + \dt J tdt 



t 2 V A 2 (A-l) 2 (A-l) 3 
,„ N d 2 A l /dA\ 2 IdA A 2 / pt , 5t 2 

(Pm) ^ = aUJ _ i^ + ^r + ^ +7A+ 4A^ 

The third equation Pni is slightly modified; the original equation can be reproduced by 
the simple change of variables (t, A) — > (t 2 ,t\). It is well known that these equations are 
characterized by the absence of "movable singularities" other than poles. 

R. Fuchs proposed two more approaches to the sixth equation P V i- One approach is 
the concept of isomonodromic deformations. In this approach, Pyi is interpreted as a dif- 
ferential equation describing isomonodromic deformations of a linear ordinary differential 
equation on the Riemann sphere. This is the origin of many subsequent researches. An- 
other approach relates Pyi to an incomplete elliptic integral. Painleve |3J] took the second 
approach, and derived a new expression of Pyi in term of the Weierstrass p-function. This 
work of Painleve is briefly reviewed in Okamoto's work on affine Weyl group symmetries 
of Pvi §■ 

Manin || revived the almost forgotten work of Fuchs and Painleve after nearly ninety 
years. Manin's remarkable idea is to use the elliptic modulus r, rather than t, as an 
independent variable. The outcome is a Hamiltonian system with a Hamiltonian of the 
normal form 7i = p 2 /2 + V(q), where the potential is a linear combination of the Weier- 
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strass p-function and its shift by three three half periods. This is a non-autonomous 
system, because the Hamiltonian depends on the "time" r through the r-dependence of 
the jp-function. 

Levin and Olshanetsky [7| pointed out that Manin's equation resembles the so called 
Calogero-Moser systems, i.e., the various extensions || of the integrable many-body sys- 
tems first discovered by Calogero M. More precisely, the Hamiltonian 7i is identical 



to a special case (the rank-one elliptic model) of Inozemtsev's extensions |fT0| , |TT|| of 
the Calogero-Moser systems. Levin and Olshanetsky called this relation the "Painleve- 
Calogero correspondence" . 

One will naturally ask if this correspondence can be extended to the other Painleve 
equations. Manin himself raised this problem in his paper. Olshanetsky jEJ conjectured 
that a degenerate version of Inozemtsev's elliptic model will emerge therein. 

This paper aims to answer this question affirmatively A guiding principle is the 
degeneration relation of the six Painleve equations ||13|| . This relation can be schematically 
expressed as follows: 



Pvi - 


- Pv - 


- Piv 






1 


1 






Pin - 


-> Pll - 


Pi 



This diagram means, for instance, that Py can be derived from Pyi by a degeneration 
process, which amounts to confluence of singular points of the aforementioned linear 
ordinary differential equation in the isomonodromic approach. We shall trace this process 
carefully on the "Calogero side" , and find a Py- version of Manin's equation. In principle, 
one can thus find an analogue of Manin's equation for all the six Painleve equations 
(though, actually, one can resort to a more direct approach that bypasses the complicated 
degeneration process). 

Remarkably (or rather naturally?), all the six equations on the Calogero side turn out 
to become a (non-autonomous) Hamiltonian system with a Hamiltonian of the normal 
form 7i = p 2 /2 + V(q). Furthermore, the Hamiltonians on the Calogero side of Py and 
Piv coincide with the Hamiltonians of the (rank one) hyperbolic and rational models 
in Inozemtsev's classification (which were discovered by Levi and Wojciechowski JL4 



before Inozemtsev's work). Those corresponding to the other three Painleve equations are 



3 



not included therein, but may be thought of as a further degeneration of the hyperbolic 
and rational models. 

One can further proceed to the higher rank models, and ask if there is still a Painleve- 
Calogero correspondence. We shall show that this is also the case. The Painleve side of 
the correspondence is a kind of multi-dimensional extensions of the Painleve equations. 
They are obviously different from another multi-dimensional extension called the "Gamier 
systems" ||13|| . For this reason, we call our mult i- dimensional extension a multi- component 



version of the Painleve equations. 

This paper is organized as follows. Section 2 is a brief review of the work of Fuchs, 
Painleve and Manin. Section 3 deals with Py, Prv and Pm- The degeneration process is 
discussed in detail for the case of Py. The direct approach is illustrated for the case of Piv 
and Pm- Section 4 shows a reformulation of the foregoing calculations in the Hamiltonian 
formalism. The status of Pn and Pj is also clarified therein. Section 5 is devoted to 
the higher rank Inozemtsev Hamiltonians and the multi-component Painleve equations. 
Section 6 is for concluding remarks. Part of technical details are gathered in Appendices. 



II Painleve- Calogero Correspondence for P 



VI 



We here briefly review the work of Fuchs, Painleve and Manin. 
Fuchs rewrites Pyi into the following form: 

" A dz 



t(l-t)C t f 

J oo 



y/X(X-l)(X-t) 



y/ Z {z-l)(z-t) 

, 0* , 7(*-l) 
ot H 

A 2 (A-l) 2 
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1\ t{t-l) 



2 J (A-t) 2 

Here Ct is the linear differential operator (Picard-Fuchs operator) 



(2) 



which also appears in the Picard-Fuchs equation of complete elliptic integrals. In this 
respect, Pyi may be thought of as an inhomogeneous (and nonlinear) analogue of the 
Picard-Fuchs equation. 

Painleve and Manin make use of a parametrization of the elliptic curve 



z{z- l){z-t) 



(3) 



by the Weierstrass p-function. Let p{u) be the p-function with primitive periods 1 and 
t: 



p(u) = p(u I 1,t) = \ + V (7 — ; — — Z — ^2)' 

u 2 \ (u + m + nrY (m + nrY I 



(4) 



(m,n)^(0,0) 

The parametrization is now given by 



_ g(u) - ei p'H 

ea-ei ' y " 2(e 2 - ei ) s / 2 ' 1 j 

where e„ = p(u> n ), n = 1,2,3 are the values of p{u) at the three half period points 

wi = 1/2, w 2 = -(1 + r)/2, cu 3 = r/2. 

Manin's excellent idea is to do a simultaneous change of the dependent variable A — > g 

by 

A = (6) 
e 2 - ei 

and the independent vrariable £ — > r by 

* = (7) 
e 2 - ei 

Manin presents the beautiful formula 

— = - (8) 

dt *(t-l)(e 2 -ei)' 1 ; 

for the Jacobian of the latter, which plays a key role in his calculations. Pyi is thereby 

mapped to the equation 

(2m) 2 ^ = J2^nP(q + u; n ), (9) 

ra=0 

where the parameters on the right hand side are connected with the parameters of Pyi as 
a — a, a\ — —j3, a 2 = 7, a 3 = — 5 + 1/2. This equation is equivalent to the Hamiltonian 
system 

dq &H dp &H , , 

dr op dr oq 



with the Hamiltonian 



2 3 

H = \-Y, a Mq + w n ). (11) 



2 

n=0 
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Ill Correspondence for Py, Prv an d P 



in 



III.l Degeneration of Pyi to P 



V 

The degeneration of Pyi to Py is achieved by rescaling the time variable and the param- 
eters as 

t = l + et, a = a, (3 = (3, 7 = 2 5 = - (12) 

e e l e A 

and letting e — ► while leaving a, . . . ,7 and t finite [ |I3|| . 

The building blocks of Fuchs' equation ([|) turn out to survive this scaling limit as 

follows: 



~o d 2 ~d ( ~d^ 



1. The Picard- Fuchs operator: 

til - t)C t — * t 2 ^- + t^ = [t^ 
y 1 dt 2 dt V dt 

2. The sum a + • • • of four terms on the right hand side: 

, ft , l(t-l) , f x 1\ t(t-l) ft Jt 2 (X + l) 

« + TT7 + 7^ TTn + [0 — - ] 7^ > OL + — + 



A 2 (A - l) 2 V 2/ (A-t) 2 A 2 (A -l) 2 (A-l) 3 ' 

3. The quare root on the right hand side: 

y/X(X-l)(X-t) — > >/A(A-l). 

4. The incomplete elliptic integral: 

dz f x dz 



y/z(z-l)(z-t) ioo^-1)' 

In particular, the degeneration of Pyi to Py is associated with the degeneration of the 
elliptic curve to a rational curve, 

y 2 = Z ( Z -l)( Z -t) — y 2 = ^ - l) 2 , (13) 

or, equivalently, the degeneration of the torus C/(Z + rZ) to the cylinder C/Z. 

Thus, rewriting a, (3, 7, 5 and t to a, /3, 7, 5 and £, we obtain the following equation as 
a Py-version of Fuchs' equation: 

d\ 2 f x dz r-^ 1X / /3 7t 5t 2 (X + l)\ 



III. 2 Analogue of Manin's equation for Py 

As an counterpart of the q- variable for Pyi, we now consider 

f x dz 

q= L7^ry (15) 

If one prefers to being more faithful to Manin's parametrization, one should rather define 
q as 

dz 

q 



2m Jz{z - 1)' 

because 2(e2 — ei) 1 / 2 — > 2%i as Imr — > +00 (see Appendix B). Since there is no substantial 
difference, let us take the first definition that is slightly simpler for calculations. 
Let us rewrite ( |14D in terms of q. The integral can be readily calculated as 

«-*($tt)- (16) 

so that the inverse relation can be written 

V\ = -coth(g/2). (17) 
Terms on the right hand side of ( |Hp can be calculated as follows: 

,f\(\ i\ 1 - sfnh(g/2) 
VA(A_1) P " ~cosh 3 (g/2)' 

v / A(A-l) ( - A _ ^ 2 = -^sinh(g), 

The differential equation for q eventually takes the form 

d\ 2 dV{q) 



where 



^ dt J dq 



sinh 2 (g/2) ~ cosh 2 (g/2) ~2~ 



= - — - " + i- cosh(g) + ^- cosh(2g). (19) 



This gives a Py-version of Manin's equation. Note that this equation can be readily 
converted to a Hamiltonian system with the Hamiltonian TC = p 2 /2 + V(q). 

Remark. 

A very similar change of dependent variable for Py is discussed in the book of Iwasaki 
et al. Il5| . 



III. 3 Idea of direct approach 



Although the degeneration process can be continued to the other Painleve equations, 
we now present a more direct approach. Note that the integrand is connected with 
the coefficient of (dX/dt) 2 in the original Painleve equation by the following very simple 
relation: 



1 



y/z(z-l)(z-t) 
1 



exp 



cxp 



1/1 I 

2 \z 
1 

2z + z - 1 



1 



z — 1 z — t 

dz 



dz 



If this is a correct prescription, one will be able to define the g-variable for Pm and Pn 
directly without the cumbersome degeneration process. This is indeed the case, as we 
shall show below. 



III. 4 g-variable for Pry 

Since the expected integrand is given by 



ex P - / 7T 



dz\ 1 



2z J ^fz 



we define 



(20) 



z 



(21) 



This can be solved for A as 



(22) 
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Honest calculations show that all derivative terms of Prv can be absorbed by the second 
derivative of q: 

d 2 q 1 d 2 X _ 1 fd\\ 2 

~d& ~j\d& 2\VX \dt) 

= 7= x {l x3 + m2 + 2{t2 - a)x + \)- (23) 

Substituting A = (q/2) 2 gives the second order differential equation 

d 2 q dV{q) 



dt 2 ' dq 

with the potential 



(24) 



^)-5(I)- 2f (l)-^-»)(f)^(I) <-) 



III. 5 g-variable for Pni 

The integrand is expected to be given by 



exp(-/f) 4 (26) 

We consider 

C x dz 

q = J - = logA (27) 

and its inversion 

A = e q . (28) 

All derivatives terms of Pni are now absorbed by the second derivative of q with respect 
to logt: 

d\ 2 t 2 d 2 X t dX t 2 fdX x 1 



t dt) 9 X dt 2 + Xdt X 2 \dt 



aX fit 7A 2 St 2 , , 

TT + TV + V + 7T2 • 29 
4 4A 4 4A 2 



Substituting A = e q gives the second order equation 



with the potential 



V(q) = -\e* + f e- - |e 2 « + ^ . (31) 
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III. 6 Summary 

Let us summarize the results of this section. 

Theorem 1 The foregoing change of variable A — > q maps Py ; Prv an d Pm to a second 
order differential equation for the new dependent variable q. These equations are equiv- 
alent to a non- autonomous Hamiltonian system with a Hamiltonian of the normal form 
H = p 2 /2 + V(q): 

(Pv) The Hamiltonian system takes the form 

dq dH dp dH 
t dt = W t tt = ~~dq~ (32) 

with the Hamiltonian 

p 2 a (3 yt . , . 5t 2 , . 

H = 7 ~ ^75) " ^J2) + T cosh ("' + T aM M- (33) 

(Piv) The Hamiltonian system takes the form 

dq dH dp OH 

dt dp ' dt dq 
with the Hamiltonian 

H= p 4-l(iY-2t ay - 2 { t 2 - a) a)\p ay 2 . m 



(34) 



2 \2J \2J v ' \2J \2 



(Pm) The Hamiltonian system takes the form 

dq dTi d% 
dt dp ' dt dq 



4 = ^ 4 = -^ (36) 



with the Hamiltonian 

" = T-? + f e ~*-I e2 ' + f e " 2 '- (37) 

Remark. 

1. The Hamiltonians for Py and Prv coincide with those of the hyperbolic and rational 



models of Inozemtsev ||10|| , Levi and Wojciechowski The Hamiltonian for Pm 



has no counterpart in their work, but nowadays can be found in the literature JT(| . 

2. The foregoing construction of the q- variable does not literally work for Pn and Pi, 
because there is no (dX/dt) 2 term. The status of these equations will be clarified in 
the next section from a different point of view. 
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IV Hamiltonian formalism of correspondence 



IV. 1 Hamiltonians of Painleve equations 

All the six Painleve equations are known to be expressed in the Hamiltonian form 

dX _dH dX _ dH 
dt dfj, ' dt dX 

with a suitable choice of the canonical conjugate variable /i and the Hamiltonian H |18| . 

This expression is by no means unique; we here consider the following Hamiltonians |13|] . 

These Hamiltonians are referred to as the "polynomial Hamiltonians" because they are 

polynomials in A and /z: 



(Pvi) H 

(Pv) H 

(P IV ) H = 2X 

(Pin) H = ^ 

(Pn) H 



A(A - 1)(A - t) 



t(t-l 
A(A-l) 2 



2 / K . K l 

A 1 ~ T + 



9-1 



K 



+ 



A A- 1 X-t 



A A — 1 (A-l) 



A(A-l) 

K 



A(A-l) 



X K \ 6»oo 

— h t + — u H 

2 X J 2 



7 _o ^Vot 
X X 



V - \Voc + - + 



^00(^0 + goo) 

2A 



2 



n + - I A. 



(Pi) # = ^--2A 3 -tA. 

Here Ko, «i, etc. are constants that are connected with the parameters a, j3, 7, 6 of the 
Painleve equations by simple algebraic relations: 



(Pvi) a 



(Pv) « 



{K + Ki + e- r 



r«o + Si 



(Prv) a = 2^ - + 1, = -2/cg. 

(P m ) a = /3 = 4r] (9 + 1), 7 = 4t&, 5 = -Arjl 

IV. 2 How to find canonical transformations 

The goal of this section is to show that the Painleve- Calogero correspondence is, in fact, 
a (time-dependent) canonical transformation of two Hamiltonian systems. By this, we 
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mean that the functional relation between A and q can be extended to (A, /i) and (q,p) so 
as to satisfy the equation 

p,dX — Hdt = constant • (pdq — TidT) + exact form. (38) 

with a suitably redefined time variable T (such as the logarithmic time logt in Py and 
Pin)- The constant factor on the right hand side is inserted simply for convenience; if 
necessary, one can normalize the constant to 1 by suitably rescaling p, q, 7i and T. For 
this reason, wel call this type of coordinate transformation a "canonical" transformation 
even if the constant factor is not equal to 1. 

Let us illustrate, in the case of Pyi, how to find such a canonical transformation. 
Suppose that A and /x be a solution of Pyi in the aforementioned Hamiltonian formalism, 
and that q be a corresponding solution of Manin's equation. The canonical equation for 
A takes the form 

dX _ A(A-l)(A-t) / _ ko _ ki _ 0-l \ 
~dt ~ t(t - 1) V ^ _ T ~ A - 1 ~ X-t) ' 

This equations can be solve for \i: 

t(t - 1) dX 1 /«o Ki 9-l\ 
^ = 2A(A-l)(A-t) ~dt + 2 \~X + A- 1 + X-t ) ' 

Our task is to rewrite the right hand side in terms of p and q. We first consider dX/dt. 
Differentiating @ against t gives 

= ( p'(q) d( l + ^ / A dr 
dt \€2 — ei dr T ) dt } 

where we have introduced the functions 

/W = » M" = — ^ — • ( 39 ) 

e 2 - ei dr 

The derivative dq/dr can be read off from the canonical equation for q: 

dq 1 dH p 

dr 2ni dp 2m 

As for the Jacobian dr/dt, Manin's formula (§) is available. One can thus express dX/dt 
as a function of p, q and r. The other part of the foregoing expression of fi contains A 
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only, which can be readily converted to a function of q and r by We thus obtain the 
following expression of /x: 

e 2 - ei 27u(e 2 - ei) 2 
^ = ,/ n P+ fr{q) 

eo — ei / /tn — 1 \ . , 

I 1 1 . (40) 

2 \p(q)-ei p(q)~e 2 p(q)-e 3 J 

We now move the point of view, and think of (^) and ( [Hp as defining a coordinate 
transformation (A, fi) — > (q,p). This gives a canonical transformation that we have sought 
for: 

Theorem 2 and gQ) defi ne a canonical transformation that connects the Hamilto- 
nian form of Pyi and Manin's Hamiltonian system. The canonical coordinates and the 
Hamiltonians of the two systems obey the equation 

dr 

fidX — Hdt = pdq — + exact form. (41) 

TV. 3 Proof of Theorem @ 

Total differential of (|6|) gives 

d\ = dq + f T (q)dr, 

e 2 - ei 

so that fidX can be expressed as 

V p (?) p (q) J \e 2 - ei 

+ 2(T + A^T + X^r A 



where 



pdg+(A) + (B) + (C), 



A = \ - /r(g)dg, 
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As we shall prove in Appendix A, (A) can be further rewritten 



(A) 



p{q + ^3) 

Am' 



7T 



fAq) 

/'(?) 



dr + exact form, 



where f'(u) denotes the w-derivative of f(u): 

df(u) p'{u) 



For (B) and (C), we have 



(B) 

(C) 

Thus we find that 



/'(«) 



p + 2m 



du 



e 2 - e x 



(42) 



(43) 



fAq) 

/'(?) 



dr. 



f'(q) 

2{\-t) dt + \ (K ° bg X + Kl bg(A - X ) + - X ) lQ g( A " *)) 
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2(X-t) 



dt + exact form. 



~ dr 

/j,d\ — Hdt = pdq — 7i h exact form, 

2m 



where 



dt , 

H = 2m — ( H 

dr 



9-1 

2{X-t) 



2vrz 



p(<? + ^s) . f T {q) . ■ 



47TZ 



fAq) 

f'(q) 



(44) 



(45) 



Our task is to prove that the transformed Hamiltonian 7i coincides, modulo irrelevant 
terms, with the Hamiltonian of Manin's equation. Here "irrelevant" means that the term 
is a function oft only. Such a "non-dynamical" term can be absorbed by the "exact form" 
part of the foregoing relation of 1-forms, thereby being negligible. 

Let us evaluate the contribution of 2m\dt/ dr)H . By Manin's formula (§) of dr/dt, 
and also by the identity 



A(A - 1)(A - t) 

we can rewrite 2m\dtj 'dr)H as follows: 

P'(q) 2 



4(e 2 - ei )3' 



n dt TT 

2m— H 

dr 



2(e 2 - ei ) 2 

P\q? 
2(e 2 - ei ) 2 

p\q? 



A A — 1 ' X-t 



H + 



A(A-1)J 



9-1 



X A - 1 X-t 



2(e 2 -ei 



Ki 9-1 

+ 



4 V A A- 1 X-t 



+ 



A(A-l) 
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The first term on the right hand side is equal to 



1 



- [p + 2m- 



JAq) 



/'(?) 



P 



JAq) 



/'(?) 



JAq) 



/'(?) 



by which the terms proportional to f T (q)/ f'(q) and its square in the definition of 7^ are 
cancelled out. The transformed Hamiltonian H can now be expressed as 



p'{qf 



(0-l)*(t-l)(e 2 -ei) 



2(e 2 -ei 



A-t 



2(e 2 - ei 



1 / /t 



0-1 



4 V A A - 1 A - t 



+ 



ft 



A(A-l) 



(46) 



Note that this is already of the normal form p 2 /2 + V(q) with the potential 

p'{qf (0-l)t(t-l)( ea - ei ) 



2(e 2 - ei ) 2 
2(e 2 - ei ) 2 



A-t 

ftn K\ 9 — 1 

— H 1 

A A- 1 X-t 



+ 



ft 



A(A-l) 



(47) 



What remains is to express V(q) as an explicit function of q. To this end, we substitute 
the factor p'(g) 2 /2(e 2 — e^ 2 by 2(e 2 — ei)A(A — 1)(A — i), and rewrite the main part of V(q) 
as a linear combination of A, 1/A, 1/(A — 1) and 1/(A — t). This leads to the following 
expression of V(q): 

V(q) = (e 2 -ei)A 

_«g (eg - e l )t k\ (e 2 - ei )(l-t) (0 - l) 2 + 1 (e 2 - ei )t(t - 1) 
2' A 2' A-l 2 ' X-t 

— -p(q + CO3) + function of t only. 

The final piece of the ring is the general formula 

[ej -e k ){ej -e e ) 



p(u + UJj) = 6j + 



(48) 



where (j, k, I) is a cyclic permutation of (1, 2, 3). This implies that 

(e 2 - ei)i 



A 

(e 2 -ei)(l-t) 

A - 1 
(eg - ei)t(t - 1) 
A-t 



= p(g + ^)-e 2 , 
= p(q + ^ 3 )-e 3 , 
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so that 



~ . , (« + «i + - l) 2 - 4k . , 

^(g) = g — - + 

# 2 

— ^"^(^ + ^2) - y^(g + ^s) + function of r only. (49) 

Apart from the last term which is negligible, this potential is indeed the same as Manin's 
potential V(q) (recall the algebraic relations connecting the constants kq, etc. and the 
parameters of Pyi)- This completes the proof of the theorem. Q.E.D. 



IV. 4 Canonical transformation for Py 

This heuristic method for constructing a canonical transformation can be applied to the 
other Painleve equations. We here consider the case of Py. 

Let A be a solution of Py, fi the canonical conjugate variable, and q the corresponding 



solution of ([181) . The canonical equation for A can be written 
dX _ A(A - l) 2 ( _ ko _ Jx_ + Vit 



dt t V A A - 1 (A - 1) ; 

This equation can be solved for fi as 

1 dX 1 f k 9 1 r)xt 
At = 7WT\ T^J t ~JI + n ~T + 



2A(A-1) 2 dt 2 \ X A - 1 (A-l) 2 , 
By differentiating (|i"7|) against t and using the canonical equation tdq/dt = dH/dp = p, 
we obtain the identity 

which can be used to rewrite the expression of \x as 

V 1 (kq 61 Vit \ 

We now reinterpret ([17]) and (|5CiD as defining a coordinate transformation (A,//) — > 
{q,p)- This indeed turns out to give a canonical transformation that we have sought for: 

Theorem 3 ([77) and ( |5^ j define a canonical transformation that connects Py and the 
Py -version of Manin's Hamiltonian system. The canonical coordinates and the Hamilto- 
nians of the two systems obey the equation 

jidX — Hdt = - (pdq — TC-] + exact form. (51) 
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Proof. Since dX and dq are connected by the relation 



d\ = VX(X - l)dq, 



fid\ can be expressed as 
fidX ■ 



dX 



l pdq -2(X-l) 



lh dt+ l -d (k log A + 6 1 log(A - 1) + 



so that 



fidX — Hdt = - (pdq — H— ] exad form. 



where 



H = 2Ht + 



Vit 

A- 1' 



We can rewrite 7i to a normal form as 
ft = 2A(A-1) 2 

\2 



1 (kq , 01 



01* 



2 VA A-l (A-l) 1 



+2A(A- 1)' 



0i 



4 V A A-l (A-l) 2 



+ 



A(A-l) 



+ 



Vit 
A-l 



where 



v (q) = _^1>!^ + ». 



A A-l (A-l) 2 



1 



+ 2/c(A-l) + 



A-l 



Ko Q\ _ \ ____ 

- + - + K& 2 ^ sin h 2 (g/2) + Ycosh 2 (g/2) 

?7i(0i + l)i , , . rih 2 . . . 
H cosh(g) cosh(2g) + function of t only. 



(52) 



(53) 



(54) 



(55) 



Apart from the last negligible term, this coincides with the potential V(q) in the statement 
of the theorem. Q.E.D. 
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IV. 5 Canonical transformation for Pry 

We now consider the case of Prv- 

Let A be a solution of Prv, \i the canonical conjugate variable, and q the corresponding 



solution of (|24|). The canonical equation for A can be written 

dX 



4A/x- (X 2 + 2tX + 2n ), 

(J/L 



which can be solved for // as 



1 dX 1 A „ 2k 
/i = — — + - A + 2t + 



4A dt 4 V A 
By (p2|) and the canonical equation dq/dt = dH/dp = p, we have the identity 



^^P+l(x + 2t+^-). (56) 



so that 



4\/A 4 V A 

Theorem 4 and ( \5(\ ) define a canonical transformation that connects Piv and the 
Pjy-version of Manin's Hamiltonian system. The canonical coordinates and Hamiltonians 
of the two systems obey the equation 

jidX — Hdt = -(pdq — Tidt) + exact form. (57) 

Proof. Since dX and dq are connected by the relation 

dX = \Z~Xdq, 

lidX can be expressed as 

^dX = \pdq + - ( X + 2t + ^ ] dX 



so that 



4 r ' 4 \ A 

111 fx 2 

-pdq - -Xdt + -d ( — + 2t X + 2k log A 



fidX — Hdt = -(pdq — Hdt) + exact form, (5£ 
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where 



H = AH + 2A. 



(59) 



We can rewrite the transformed Hamiltonian TC to a normal form as 



n 



IX 



1 f X K 

"-2 2 +(+ T 



+ 8A 



1 fx 



4 V 2 



«o 



X 



+ 2A 



where 



V(q) 



-2A 
1 



A 

- +t + -T- 

2 A 



+ AOooX + 2A 
A 3 - 2iA 2 - 2(t 2 + k - 20OO - 1)A - 2k 2 A- 1 



+function of i only. 



(60) 



(61) 



Substituting A = (q/2) 2 gives the potential V(q) modulo an irrelevant term. Q.E.D. 



IV. 6 Canonical transformations for Pni 

The situation of Pni is somewhat similar to Py- 

Let A, again, be a solution of Pm, A the canonical conjuage variable, and q be the 
corresponding solution of (|30|). The canonical equation for A takes the form 

dx a 2 / e Vo t 

which can be solved for // as 



t dX 1 / #o Vot 



~ 2A 2 dt + 2 y 00 + A A 2 
By differentiating (p8|) and using the canonical equation tdq/dt = dTi/dp = p, the t- 
derivative of A can be written 

dX , 



so that we obtain 



^ = 2A 2r 00 a"^'' (62) 



This relation, again, can be used to define a canonical transformation: 
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Theorem 5 ( fjffij and define a canonical transformation that connects Phi and £ne 
Pni-'werszon of Manin's Hamiltonian system. The canonical coordinates and the Hamilto- 
nians of the two systems obey the equation 



jid\ — Hdt 



pdq — H 



.dt 
T 



exact form. 



Proof. Since dX and dq are connected by the relation 

dX = Xdq, 

fidX can be written 



fidX 



1 7 1 
2 pdq + 2 



do _Vot 
X X 2 



dX 



A 



so that 



where 



/idA — Hdt 



pdq — 7Y— 



exact form, 



n = 2Ht + ^4- 

X 

We can convert the transformed Hamiltonian 7i to a normal form as 

2 



2A 2 



+2A^ 



2 I X 



^ + A ~ 



' A 2 
A 2 



+ 



?7oo(flo + goo) 

2A 



+ 



where 



?7oo + 



A 



Vot 
A 2 



+ «oo(^0 + ^oo)A + 



A 



"oo^ooe" + r^o + l)*e"« - ^e 2 * - ^e"^ 
+function of t only. 



(63) 



(64) 



(65) 



(66) 



(67) 



Thus, apart from the last irrelevant term, V(q) coincides with the potential V(q) in the 
statement of the theorem. Q.E.D. 
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IV.7 Status of P n and Pi 

Let us turn to P n and Pi. The Hamiltonian of Pi is already of the normal form H = 
\ + V(q) with A = q, \x = p and H — H. Although this is not the case for P n , one can 
directly find a canonical transformation that converts the Hamiltonian H to a normal 
form: 

Theorem 6 A Pu-version of Manin's Hamiltonian system is defined by the Hamiltonian 

2 1 / j \ 2 

V 1 ( 2 , tX 



H = -2~2\ q + V <68) 
This system is connected with Pn by the canonical transformation 

X = q, f i=p + X 2 + t -. (69) 

The canonical coordinates and the Hamiltonians of the two systems obey the equation 

fid\ — Hdt = pdq — Hdt + exact form. (70) 

Proof. The foregoing relation between (A,//) and (q,p) implies that 

, A„ A , /A 3 tX\ 

IxdX = pdq + I A + - J dX = pdq — —dt + d I — + — I , 

so that 

jidX — Hdt = pdq — Hdt + exact form, (71) 



where 



H = B + \ 

2 i / . \ 2 



2 1 / j \ 2 



1 tY ( \\ . A 



^ A + - - a + - A + - 



2 V 2 / V 2 / 2 

= t-H ,2+ sJ (72) 

This is nothing but the Hamiltonian in the statement of the theorem. Q.E.D. 
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V Multi-component Painleve equations 



V.l Inozemtsev Hamiltonians of higher rank 

The rank i version of Inozemtsev's Hamiltonians have i coordinates q±, . . . ,qe and canon- 
ical conjugate momenta Pi, ■ ■ ■ ,pe- The Hamiltonians of the elliptic, hyperbolic and 
rational models take the following form JTDL O, 



Elliptic model: 

w = ^2 [ ~2 + E 9 ™ p ( q j + Un "> ) + & E ~ q ^ + p ( q j + ■ 

j=l \ n=Q J j^k 

Hyperbolic model: 

1 1 



Rational model: 

~2 



sinh 2 ((g,- - q k )/2) smh 2 (( qj + q k )/2) 



n = E (f + ^ + £<6 + gk] + ghf) + g\ E 



+ 



(ij - Qk) 2 (qj + Qk) 



Here g§, g\, gi, g?> and g 4 are coupling constants. The Painleve- Calogero correspondence 
for Pjn, Pn and Pi suggests the existence of further degeneration of these models. 

The goal of this section is to extend the the Painleve-Calogero correspondence to 
these higher rank models. Since a complete exposition will become inevitably lengthy, 
we shall illustrate the elliptic and hyperbolic models in detail, leaving the other cases 
rather sketchy. The strategy is as follows: The point of departure is the Hamiltonian of 
Inozemtsev's rank t elliptic model. This gives rise to a rank I version of Manin's equation. 
Starting with this non-autonomous Hamiltonian system, we seek for an analogue of the 
degeneration process for the Painleve equations. We can thus obtain six types of non- 
autonomous Hamiltonian systems. At each stage of the degeneration process, we confirm 
that the non-autonomous Hamiltonian system on the Calogero side can be mapped, by a 
canonical transformation, to a multi-component analogue of the Painleve equation of the 
corresponding type. 
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V.2 Elliptic model and multi-component Pyi 

We now consider the non-autonomous Hamiltonian system 

_ dqj OH .dpj dH 

ar opj clt oqj 

defined by the Hamiltonian of Inozemtsev's elliptic model. This is a rank I version of 
Manin's equation. This non-autonomous system is known to describe a family of isomon- 
odromic deformations on the torus [ |17| |. 

An honest generalization of the canonical transformation for the case of I = 1 leads 
to a mult i- component version of Pyi as follows: 



Theorem 7 The time- dependent canonical transformation defined by 
A, 



ea-ci ' 
e 2 -e 1 2ni(e 2 - e{) 2 

* = pHY fM]) 

+ e2 ~ ei ( * + ^ + °~ l ). (71) 

2 \p(<lj) ~ e i p(lj) ~ e 2 p(qj) - e 3/ 



and 



t = (75) 
e 2 - ei 

maps fffflj to the Hamiltonian system 

d\j dH dfXj dH 
dt dfij ' dt d\j 

with the Hamiltonian 



(76) 



« = £ 



A,(A,-l)(A.,-«) 



^ - ( — + ~\ 7 + x 7 ) »i 



Xj Xj - 1 Xj - 1) Xj(Xj - 1) 



g\ x - - -t) + X k (X k -l)(X k -t) 

A A j + A fcJ 



ttE 



2t « ->*» 



i(Aj - A t )= 



■ (77) 



Proof. The method of proof for the case of £ = 1 can be applied to the present case as 
well, yielding the equality 

1 dr 1 

^^Pjdqj — H- — : = fijdXj — Hdt + exact form, (78) 



2wi 

3=1 3=1 
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where 
H 



-f 



t(t-l) 



^2 _ ( ^0 + ^1 + ^ — 1 

Aj Aj - 1 Aj - t 



Hi + 



A,(A, - 1) 



2 

2t(t - l)(e 2 - ei ) ^ ( Piqj ~ Qk) + P{q > + Qk) ) ■ 



(79) 



What remains is to express the "two-body potential" part in terms of Aj. To this end, let 
us recall the addition formula 

p'{uf + p\vf 



p(u -v) + p{u + v) = -2p{u) - 2p{v) , _ 
of the p-function. Applying it to the case where (u,v) = (Aj, X k ), and substituting 

p(qj) = ei + (e 2 - e^Aj, 

p(qk) = ei + (e 2 - ei)A fc , 
(e 2 - eO 3 



0) 



A,(Aj-l)(Aj-t), 



(e 2 - ei) ; 



Afc(A fe — l)(Afc — i), 



we can rewrite the two-body potential terms as 



P(0j - Qk) + p(Qj + Qk) 



-2(ei + (e 2 - ei)Aj) - 2(e 1 + (e 2 - ei)A fc ) 
(e 2 - ei ) 3 Aj(Aj - l)(Aj - t) + A fc (A fc - l)(A fe - t) 



8 (ei + (e 2 - ei)Aj - e x - (e 2 - ei)A fc ) 2 

= -4ei-2(e 2 -ei)(Aj + A fc ) 

e 2 - ei Aj(Aj - 1)(A, - t) + X k (X k - l)(X k - t) 



+ 



• (81) 



8 (Aj - \ k f 

The first term — 4ei is non-dynamical, thereby negligible (i.e., can be absorbed by the 
"exact form" part). Removing these terms from H, we obtain the Hamiltonian H. Q.E.D. 



V.3 Degeneration of elliptic model to hyperbolic model 

The degeneration of the elliptic model is achieved by letting Imr ^ +oo. Like the 
degeneration process from P V i to Py, this is a kind of scaling limit, namely, the coupling 
constants g n and the elliptic modulus r have to be suitably rescaled. To this end, we have 
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to understand the asymptotic behavior of the constants ei,e2,es and the p-function in 
the limit as Imr ^ +00. All necessary data are collected in Appendix B. For instance, 
the asymptotic expression of ei,e2 and imply that 



t = 1 + — — - = l + lQn 2 e nir + 0(e 2nir ). (82) 
t2 — ei 

This is indeed consistent with the scaling rule t = 1 + et in the degeneration process of 
Pvi to P v . 

Having these data, we now rescale the coupling constants and the elliptic modulus as 

2-2 2 -2 2 92 1 93 3 9z 2-2 / OQ \ 

9o = 9 , 9i=9i, 9 2 = — + -^, 9s = -^, 9a = 9a (83) 

and 

16e™ r = et, (84) 

and consider the limit as e — > while leaving g n and t finite. Note that letting e — > 
amounts to letting Imr ^ +oo. 

The asymptotic expression of and p(u+u n ) in Appendix B show that the potential 
V(q) of the elliptic model behaves as 

V(q) = f (4^ + + ^003(2^) -^003(4^)) 
^ V sm (^j) cos 2 (7r<7j) 2 Jy 8 J 7 

+ ^ E f ^27-7- — Tyi + 



sin (tt(^ - q k )) sin (tt^ + q k )) 



+function of e and i only + 0(e). 
Thus, removing negligible terms, we obtain the following Hamiltonian in the limit: 

4^f p2 ij ~9> 2 ~9W ~9Wt ~gln 2 t 2 \ 

n — > — H £ - H — -H cos(27ro 1 - 2 cos(47ro,-) 

sin 2 (7rg,) cos 2 ^,) 2 1 ^ 8 1 *>> ) 

E (^WT + ■ y (\ w ) ■ (85) 

Van (tt(^ - g fc )) sin (tt(^ + g fc ))/ 

The asymptotic expression of t determines the equation of motion in the limit. In fact, 
since 

dr 7T iri 



dt t(i-l)(e 2 -ei) (l + et)(-et)(-7r 2 + 0(e)) 
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and 

n . d „ .dt dt d / 9 ~ ^ . 9 . \ d 
2m— = 2m — — — = \2>kH + 0(e 2 ) I — , 

we find that the equations of motion take the following form: 

dt opj dt oqj 

The final step is to rescale the variables and the Hamiltonian as 

and to rename t and H to t and 7i. Let us also define the new constants 

a = -?l, p = & 7 = --, S=& (88) 
2 2 2 2 v; 

which are to be identified with the four parameters of Py. The outcome is the non- 
autonomous Hamiltonian system 

t dqj_ = &H t dpj_ = _&H 
dt dpj ' dt dqj 



with the Hamiltonian 



H = ^—— |_^ + ^cosh(^) + ^cosh(2^ 

jr[ V 2 sinh 2 (g,/2) cosh 2 (^/2) 2 yq ' J 8 ) 



smh\( qj - q k )/2) smh 2 (( qj + q k )/2) 



This gives a rank t version of the non-autonomous Hamiltonian system on the Calogero 
side of P v . Note that the Hamiltoian is essentially the same as the Hamiltonian of 
Inozemtsev's hyperbolic model, except that the effective coupling constants are now time- 
dependent. 

Remark. The foregoing prescription of scaling limit of the coupling constants and the 
elliptic modulus is reminiscent of "renormalization" in quantum field theories. In this 
analogy, one can interpret the equations of motion of the Hamiltonian system as "renor- 
malization group equations" , in which t plays the role of a "mass scale" parameter. 
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V.4 Canonical transformation to mult i- component P 



v 



Again, an honest generalization of the canonical transformation for the case of i = 1 leads 
to a mult i- component version of Py: 

Theorem 8 The time- dependent canonical transformation defined by 

1 ( k Ox 7frt 



Xj = — coth(<7j/2) 
Pj 



2^(A i -l) 2 \\j Xj-1 (Xj-iy 
maps to the Hamiltonian system 

dXj dH djjij dH 

dt d/j,j ' dt dXj 



(91) 



(92) 



with the Hamiltonian 



" = E 

3=1 



t 



X, X, - 1 (Xj - 1) 



^^2(A i -l)(A Jfe -l)(A i + A Jfc ) 
2t 



(A, - A fc )^ 



(93) 



Proof. The method of proof for the case of £ = 1 can be used as it is. The outcome is 
the equality 



dt 



^^Pjdqj — = 2 | f'jdXj — Hdt J + exact form, 



3=1 



(94) 



where 



» = E 



d=i 



Vit 



A, Xj - 1 (Aj - l) 2 
1 



Hj + 



K 



A,(A,-1) 



2t ^ Vsinh^fe - q k )/2) ' sinh 2 ((g, + g fe )/2) 

The two-body potential part can be rewritten by use of the identity 

cosh(2w) cosh(2t>) — 1 



1 1 

+ 



sinh 2 (u — v ) sinh 2 (u + v) 



4 



(95) 



(96) 



(cosh(2u) — cosh(2u)) 2 
Substituting u = qj/2, v = q k /2, and also using the equality cosh(g,,) = (Aj + l)/(Aj — 1 
we find that 

1 1 2(A i -l)(A Jfe -l)(A i + A Jfe ) 



+ 



sinh^((g i - q k )/2) sinh'((g j + q k )/2) 
which gives the two-body potential term in H . Q.E.D. 



(Xj - x k y 



(97) 
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V.5 Other models 

The degeneration process can be further continued, and leads to four more models that 
correspond to a multi-component version of Ppy, Pm, Pn and Pi. Since the details of 
derivation are more or less parallel, we show the final results only. The Hamiltonian 
of each model, like those in the foregoing cases, becomes a sum of i copies of the one- 
component Hamiltonian and Calogero-like two-body potential terms. 

V.5.1 Rational model and multi-component P IV 

This model can be derived from the hyperbolic model by degeneration. The degeneration 
process consists of putting the variables and the parameters as 



t — 1 + 2et, qj = iri + e 1 ^ 2 ^, pj = 



Pj 



(98) 



2e l/2> 



and 




(100) 



(99) 




5 (?)'-» (?)' -*■-> (?)'+'(?) 



-2 




(101) 



The canonical transformation defined by 




(102) 



maps the foregoing non-autonomous system to the Hamiltonian system 

d\j dH djjbj dH 



(103) 



dt d/ij ' dt d\j 

with the Hamiltonian 



( 




2(A,- + A fc ) 
(A,-A fc ) 2 ' 



(104) 



28 



V.5.2 Exponential- hyperbolic model and multi-component Pni 

This model, too, can be derived from the hyperbolic model by degeneration. This degen- 
eration is achieved by the putting the variables and the parameters as 

Qj = -Qj ~ lo S |> PJ = -Pii ( 105 ) 

and 

and letting e — > 0. 

The equations of motion of this model takes the canonical form 

at opj at oq,j 



with the Hamiltonian 



'P] a Q , , Pt _ qj _ l 2qj ft* 



4 8 8 



-^4 2 E • 1 VJ^y ( 108 ) 



The canonical transformation defined by 

^ * = ^ + K<- + |-f) <"*» 

maps the foregoing onn-autonomous system to the Hamiltonian system 

d\j dH d/ij dH 
dt d/ij ' dt d\j 

with the Hamiltonian 



2 * t£ ( A i - A ^)' 



(110) 



H ~2^T ^-l^ + r- a^J^ + 2A~ + w2wx.._x./i2- (1U) 



V.5.3 Second rational model and multi-component P n 

This model can be derived from both the rational model and the exponential-hyperbolic 
model by degeneration. For the degeneration from the rational model, we write the 
variables and the parameters as 

_ -i +4 -V3 £ 4 t - q . _ i + 2 -i/3 £ 2g. 4 2 /% 

*- e > 2 ~ C 3/2 ' e V2 ^ 
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and 



and let e — > 0. The degeneration from the exponential-hyperbolic model is 
achieved by putting 

Pi 



t = 1 + 2e t, 5j = 2egj, pj 



and 



1 



_J_ ,/? - 1 + 4e3 ^ = J_ 
a - 2e 6' ^~ 2e 6 ' 7 ~4e 6 ' 4e 6 ' 

and again letting e — > 0. 

The equations of motion of this model takes the canonical form 

an 



dqj dH dpj 
dt dpj ' dt 



dq 3 



with the Hamiltonian 

i 

" = E 

j'=i 



1 / 2 * 



The canonical transformation defined by 



X j = q j , n j = Pj + \ 2 j + - 
maps the foregoing non-autonomous system to the Hamiltonian system 



dH dfij 



dt dfij ' dt 



dH 



with the Hamiltonian 



j'=i 



a + - ) Aj 



(A* -A*) 2 ' 



(113) 
similarly 

(114) 
(115) 



116) 



(117) 



(118) 



(119) 



(120) 



V.5.4 Multi-component Pj 

This model can be derived from the second rational model, and takes the same form on 
both the Painleve and Calogero sides. The degeneration process is achieved by putting 



^ -6 + e 12 t 1 + e% pj 15 

t= e io » * = — p ^ = 7' a = 4e 



(121) 
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and letting e — > 0. The equations of motion takes the canonical form 

dqj dH dpj dH 
dt dpj ' dt dqj 

with the Hamiltonian 



;i22) 



(123) 



VI Concluding remarks 



We have shown that the Painleve-Calogero correspondence persists for all the six Painleve 
equations and their multi-component generalizations. The Calogero side of this corre- 
spondence is a non-autonomous version of Inozemtsev's elliptic model and its various 
degenerations. Those for Py and Piv are a non-autonomous version of Inozemtsev's hy- 
perbolic and rational models. The others corresponding to Pni, Pn and Pi are further 
degenerations of the hyperbolic and rational models. The pattern of degeneration on the 
Calogero side repeats the degeneration diagram 

Pvi — ► Pv — > Piv 
I I 
Pin — > Pn — > Pi 

of the Painleve equations. 

This picture applies to the autonomous systems as well. Actually such degeneration 
relations in the autonomous case have been more or less well known to experts of Calogero- 



Moser systems (see the Introduction of van Diejen's paper |T6|). The autonomous systems 
are defined by a Hamiltonian of the same form with the time-dependent coupling con- 
stants being replaced by absolute constants (except for the elliptic model, in which case 
an independent time variable is introduced). Those in the position of the first row of 
the degeneration diagram are, of course, Inozemtsev's elliptic, hyperbolic and rational 
models (see Section 5). Those in the position of Pm and Pn are defined by the following 
Hamiltonians: 



Exponential-hyperbolic model: 

H = E (f + & qj + + + &~ 2qj ) + E 



sinh ((qj - q k )/2) 
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Second rational model: 



■ , X - / . - 1, Qk) 

3=1 v ' j^k ylJ 

The Hamiltonian in the position of Pi is redundant in the automonous case, because it is 
a specialization, rather than a degeneration, of the last Hamiltonian. 

Note that the Hamiltonian of the second rational model is a quartic perturbation of 
the usual {At type) rational Calogero Hamiltonian. According to recent work of Caseiro, 



Frangoise and Sasaki [BJ, such a quartic (integrable) perturbation always exists for any 
rational Calogero-Moser system. Inozemtsev's rational model, which is a sextic perturba- 
tion of the Dp type rational Calogero-Moser system, might admit a similar interpretation. 

Back to the Painleve equations, the extended Painleve-Calogero correspondence raises 
many interesting problems. A central issue will be to find an isomonodromic description 
of the multi-component Painleve equations. If such an isomonodromic description does 
exist, it should be related to a new geometric structure. 
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A Proof of ( gg ) 



Let us introduce the two auxiliary functions 



( \ llM hi \ + ^0 (a i \ 

9{u) = 7^y h{u) = W^Y ( } 

associated with the function 

f(u) = (A.2) 
and the standard elliptic theta function 

oo 

d{u) = exp(7iirn 2 + 2mnu). (A. 3) 
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Lemma 1 g(u) is a meromorphic function on the u-plane with additive quasi-periodicity 

g(u+l)=g(u), g(u + r)=g(u)-l. (A A) 

All poles are of the first order and contained in the lattice u 3 + Z + rZ. Furthermore, g(u) 
has zeros at u = and u — uj\. 

Proof. Since f(u) is a doubly periodic function with primitive periods 1 and r, f'(u) and 
f T (u) transform as 

f'(u+l) = f'( U ), /'(« + t) =/'(«), 

f T (u + 1) = f T (u), f T (u + r) = f T (u) - f{u) 

under the shift by 1 and r. This implies the additive quasi-periodicity of g(u). Further- 
more, by the construction, g{u) is a meromorphic function on the w-plane, and all possible 
poles are of the first order and located at the points of oo^ + Z + rZ. Let us examine the 
behavior of g(u) at the representative points u = w ,wi, —002,^3- 

• As u — > u = 0, 



/(«)=/ — ^ + 0(1) 
(e 2 - ei)it 2 



thereby 



/'(«) = - 7 ^ + 0(1), f T (u) = - e2 ' T e ,Y 2 +0(l), 

(e 2 - eijw 3 (e 2 - ei) 2 u 2 

so that has rather a zero at u — 0: 

</(«)= O(u). (A.5) 

• As u — > u>i — |, 

/(«) = — - — - ei + - wi) + 0((u - uji) 2 )] 

e 2 — ei V / 

= O((u-o*) 2 ), 

thereby 

/'(«) = 0(u - cu 1 ), f T (u) = 0((u - urf), 
so that g{u) has another zero at u — u>i. 

g(u) = 0(u-iu 1 ). (A.6) 
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As u — > -cu 2 = \ + §, 



/(«) = fp(-cJ 2 )-e 1 + p'(-cJ 2 )(M + cu 2 )+0((M + ^ 2 ) 2 ) N ) 

e 2 — ei v / 

= 0((« + cu 2 ) 2 ), 



thereby 

/'(«) = 0(« + U 2 ), f T (u) = 0(U + U 2 ), 

so that g(u) behaves as 

g{u) = 0(1). (A.7) 



• As u — > C^3 - ' 



2 ' 



/(«) = p(cJ 3 )-e 1 + p / (w 3 )(M-^3)+0((M-W3) ) 

e 2 — e\ \ 

= t + 0((u-LU 3 ) 2 ), 



thereby 



/'(«) = 0(«-ws), / T («)=0(l), 
so that ^(u) turns out to have a pole of the first order at u = cj 3 : 

( ? H = 0(( M -cu 3 )- 1 ). (A.8) 

The behavior of g(u) at the other points of uo n + Z + tZ can be deduced from these results 
by the additive quasi-periodicity of g(u). Q.E.D. 



Lemma 2 h(u) is a meromorphic function on the u-plane with additive quasi-periodicity 

h(u + 1) = h(u), h(u + t) = h{u) - 2ni. (A.9) 

All poles are of the first order and contained in the lattice u 3 + Z + rZ. Furthermore, h(u) 
has zeros at u = and u = uj\ . 



Proof. Let us recall the fundamental properties of 
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• $(u) is an entire function on the w-plane with zeros of the first order at the lattice 
points u 2 + m + nr (to, nGZ). 

• is quasi-periodic, 

0(u + 1) = 0(u), 0(u + r) = e- 7 ™"- 2 ™^). 

• and d(u + 1/2) are even under the reflection u — > — w. 

All the properties of h(u) in the statement of the lemma are an immediate consequence 
of these properties of Q.E.D. 

Lemma 3 The function f(u) satisfies the equation 

J T (u) _ &(u + u 1 ) 
m f'{u) diu + uj,)' 

where the prime stands for d/du. 

Proof. The foregoing properties of g{u) and h(u) imply the following: 

• 2mg(u) — h(u) is a doubly periodic meromorphic function with fundamental period 
1 and r. 

• All poles of 2nig(u) — h(u) are of the first order and contained in the lattice oj^ + 
Z + rZ. 

• 2irig(u) — h{u) has zeros at u = and u = uj\. 

The first two properties imply that 2nig(u) — h(u) is a constant. By the last one, this 
constant has to be zero. We thus find that 2mg{u) — h(u) = 0. Q.E.D. 

Lemma 4 -&{u) satisfies the equation 

(\og-&{u + = —p(u + Us) + function of r only. (A.ll) 

Proof. The aforementioned complex analytic properties of $(u) imply the following: 
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(A.10) 



(jog$(tt + 0Ji)\ is a doubly periodic meromorphic function with primitive period 
1 and r. 

All poles of this meromorphic function are contained in the lattice uj^ + Z + rZ. 



As u — ► — co>3, this function behaves as 

log + c^) 



(li + LU 3 ) 2 



+ 0(1). 



The function — p(u + CJ3), too, has these properties. Accordingly, their difference is a 
constant function on the w-plane, namely, a function of r only. Q.E.D. 



We now return to the proof of (S2f) . By the third lemma, we have the identity 



27r .fr(v± du _ + uj x ) ^ _ d${u + wi) <9tf(u + wi)/9r 



On the other hand, the well known "heat equation' : 



t?(u + o;i) 



(A.12) 



Ani- 

implies that 

&&{u + u x )/dT 1 ${u + oj x )" 



dr 



1 

47TZ 



#(u) 



(A.13) 



log$(tt + Ct»i] 



By the third and forth lemmas, the last line can be rewritten 

2 



" (d'{u + uji) 



0(lt + Wl) 



1 

47TZ 



-p( w + u, 3 )+ (2^^ 



/'(«) 



+ function of r only, 



so that 



9 -^ U h 1 

27T2 cfat = — 

/'(«) 47TZ 



p{U +U 3 ) - 27TZ — 



/'(«) 



(ir + exact form. 



(A. 14) 



Substituting w = q gives (|42|) 
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B Asymptotics of elliptic functions 



The asymptotic behavior of the p-function p(u), the shifted p-functions p(u + ujk) and 
the constants eu = p(uk), in the limit as Imr — > +00, can be deduced from the well 
known formula 



00 2 2 00 n 2 

PN = ^ sm\n(u + nr)) ~ T " sin 2 (vmr)- ' ' 1 



n=— 00 



Let us first consider the asymptotic behavior of p(u) itself. The constant (n = 0) term 
in the first sum is of order 1 and the n-th term is of order e 2nmT , Similarly, the n-th term 
in the second sum is of order e 2nmT . Therefore 

p(u) = -^—-^ + 0(e 2 ™ ) - (B.2) 
sin (iru) 3 

A similar estimate leads to the following asymptotic expression for the shifted p- 
functions: 

cos 2 (iiu) 3 

7T 2 

p(u + u 2 ) = -— + 8ir 2 cos(2iru)e mT + 0(e 2mT ), 

O 

7T 2 

p(u + uj 3 ) = -—-8n 2 cos(2iru)e 2nlT + 0(e 2mr ). (B.3) 
o 

In fact, the degeneration process of the elliptic model requires us to know the asymptotic 
expression of p(u +u 2 ) + p(u+ u 3 ) to the order e 2mr . This can be achieved by the following 
calculations: 

p(u + u 2 ) + p(u + u 3 ) 

7T 2 27T 2 _ y _ ^ 



cos 2 (w + I + nr) sin (u + | + nr) 3 ^-^ sin (7rnr) 
- 32tt 2 cos{2TTu)e 2mT + 16vr 2 e 2 ^ r + 0(e 37nr ). (B.4) 



3 

We now consider the constants e^. For instance, e 1 can be written 

00 2 2 00 r> 2 

x - 7T 7T x - 27T 

^ cos 2 (7rnr) 3 ^—f sin 2 (7rnr) 



?i=— 00 



= ^ 2 +i;— ^-z-^-^ (b.5) 

3 ' cos^(7rnTl f sin rar 

n=l v ' n=l v ' 
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The constant 27r 2 /3 becomes the leading term; the leading (n = 1) terms of the last two 
series give the next-leading term of the order e 2ntT . e 2 and can be similarly analyzed. 
Thus the following asymptotic formulas are obtained: 

ei = ^ + 167r 2 e 2 ^ r + 0(e 47rir ), 

7T 2 

e 2 = -— + 8n 2 e niT + 0(e 2mr ) 

e 3 = -87rV ir + 0(e 27rlr ). (B.6) 

o 

In particular, e 2 — ei — > — 7r 2 , as expected. 
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